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I. INTRODUCTION 

Limited random walks, with sudden termination of a 
trajectory are frequently analyzed in descriptions of mo¬ 
tion in porous media, biological tissues, composite mate¬ 
rials, and dynamic networks and of extreme, catastrophic 
events like gambler’s ruin, chemical reactions, and species 
extinction m- Quite often, however, the absorption 
events and disappearance of trajectories are followed by 
resets or restart activities of the system, e.g. the reloca¬ 
tion of searching paths in animal foraging, the seeking for 
target location by repair proteins or the returning to the 
initial position after an unsuccessful search of the address 
by an individual lost in a vast city [SHE]. 

A random walk with restart is also known as a graph 
mining technique widely used in the machine learning 
community for page-ranking or web search models and 
cryptology ima In this approach the frequency of visits 
paid to a given node can be analyzed as a random walk on 
a graph. It is described as an ordered sequence of visits to 
vertices with a source (initial) vertex probability po- For 
Markov chain models of transitions between subsequent 
locations on a graph described by a matrix 11, reset events 
inject additional randomness to the walk pt+i = (1 — 
c)Ilpi -\- cSi with c being the probability of resetting per 
step and Si representing an arbitrary probability vector 
added at resetting. 

Many intriguing facets of the process in which a Brow¬ 
nian particle is stochastically reset to its initial position 
with a constant rate have been investigated by Evans 
and Majumdar [5]. The stationary state of such a pro¬ 
cess has been shown to be described by a non-Gaussian 
distribution which, due to a non-vanishing steady state 
current directed towards the resetting position, violates 
the detailed-balance condition. The temporal relaxation 
towards this nonequilibrium steady state has been shown 
to exhibit a dynamical transition m- Moreover, it has 


been proved [5] that there exists an optimal resetting 
rate that minimizes the average hitting time to the tar¬ 
get. Extensions to space depending rate, resetting to a 
random position with a given distribution and to a spa¬ 
tial distribution of the target have been also considered in 
Ref. |14j . Brownian diffusion in external potentials have 
been further analyzed in a recent study by Pal |15j . 

In a somewhat different context, similar random walks 
with stochastic resets have been analyzed by Durang et 
al. |16j who posed the problem of interacting particles 
subject to a stochastic return to the initial configuration 
in the coagulation-diffusion process. The particles per¬ 
form random hoppings to nearest-neighbour sites such 
that upon the encounter of two particles, the arriving par¬ 
ticle disappears. The stochastic reset is described then 
by a given set of probabilities for having some consecu¬ 
tive empty sites. A Markov monotonic continuous time 
random walk model in the presence of a drift and Poisson 
resetting events has been addressed in an elegant work by 
Montero and Villarroel m, who derived general formu¬ 
las for the survival probability and the mean exit time. 

While most of the works related to random walks with 
resets is based on continuous and discretized version of 
a Wiener process, relatively few studies have been de¬ 
voted to resetting accompanying generalized Wiener mo¬ 
tion with discountinuous Levy jumps. Levy flights and 
Levy walks m have been claimed to be observed in many 
foraging animal species [T9H28] . which has led to the¬ 
oretical analysis showing an optimality of Levy flights 
or Levy walks in different setups [IHl . The sum¬ 

mary of those results can be found in a recently pub¬ 
lished book [34] . The optimization of a mean first pas¬ 
sage time (MFPT) in a discrete time model of Levy flights 
with stochasting resetting has been addressed in Ref. [3S] , 
where it has been shown that the optimal parameters ad¬ 
mit jumps (i.e. discontinuous changes) as functions of a 
distance to the target. Hereafter, by analyzing statistics 
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of first arrival times (FAT) of the continuous time version 
of the model, we demonstrate parameter-dependent tran¬ 
sition between the optimal Gaussian and non-Gaussian 
search strategies. 

In this work we concentrate on a variant of the model, 
in which a one-dimensional jump-like searching process 
with resetting events is analyzed as a renewal Markov 
model with Levy jumps. We assume that a random 
searcher starts its motion at xq = 0 and tries to find 
the object located at some position x. The walker does 
not memorize its former locations, and the steps under¬ 
taken at any instant in time are statistically independent 
and drawn from a symmetric stable distribution with a 
stability index a S (0,2]. Furthermore, at random times 
following the Poisson point process, the searcher decides 
to instantaneously reset to the initial position. We de¬ 
rive an expression for the transition probability density 
of such process, analyze the existence and character of 
the long-time stationary distribution and discuss optimal 
conditions for the mean first arrival time (MFAT). 

The paper is organized as follows: Section II introduces 
the model and discusses the structure and stationary so¬ 
lutions of evolution equations for corresponding proba¬ 
bility distribution functions. The mean first arrival time 
in the model is introduced in Section III and its optimiza¬ 
tion is further analyzed in Section IV which presents the 
most important results of this work. We summarize the 
paper and add conclusions in Section V. 


II. TIME EVOLUTION AND TRANSITION 
PROBABILITY 

We start with an analysis of the integral equation that 
governs the evolution of the probability density function 
for the process > 0}: 

W{x,t\xo,to)dx = 

Prob{a: < X{t) <x + dx\X{tQ) = a:o} (1) 

In the course of time IT(a:, tjxo, to) is subject to possi¬ 
ble reset events to a: = 0 or jumps (Levy flights). Re¬ 
sets are independent from flights and occurring in time 
according to Poisson statistics with an average expecta¬ 
tion time for the occurrence of the event given by r~^. 
Note that for the purpose of analysis, we have untied 
the initial and resetting positions. We denote the for¬ 
mer as Xq and keep the latter at the origin. The over¬ 
all process is time homogeneous, i.e. IT(a:, t|a:o, to) = 
W(x,t — tol^^OjO) = W{x,t — to|a^o), so that the prop¬ 
agator satisfies the equation (for the derivation and a 
detailed discussion see Ref.pH]!: 

W{x,t\xo) = e~'^^Wo{x,t\xo) + 

+ f dTe~'^'^rWo(x,T\0) ( 2 ) 

Jo 


without resetting events, whereas the second term de¬ 
scribes the evolution after the last reset. The function 
Wo(a:, rjxo, to) denotes the probability density function 
(PDF) of the process when the resetting mechanism is 
switched off. In this case the random walk propagator 
fulfills equation 


Wo(x,t|a;o) = 6{x - xq) 


1 — / $(r)dr) 

Jo 

pt P+OO 

+ / $(t — t') / p{x — x')Wo{x',t'\xo)dx'dt', 
Jo J — oo 


( 3 ) 


where <I>(t) is the waiting time PDF, independent of the 
jump-length PDF p{x — x'). In the Fourier-Laplace space 


W {k, s) = F[C\W {x,t);t ^ s]] X ^ k], 


the integral Eq.([^ takes the form of 


Wo{k,s\xo) 


1 - d>(s) 1 

s l — ^{s)p{k)’ 


( 4 ) 


where C[f{t)] = exp(—st)/(t)dt. We further assume 

that $(t) has a well defined mean value. To = t^{t)dt, 
and p{x) is the PDF of the Levy stable form, so that its 
characteristic function reads 


X[p{x)\ = exp [-CT“|fc|“], (5) 

with the stability index 0 < a < 2. The resulting process 
is Markovian, with the variance diverging for a < 2 and 
fractional moments m scaling like: 

(|x(t)r) cx (Dt)^/“, (6) 

where D = a°‘lTQ. The asymptotic behavior of 
Wo(fc, sjxo, 0) can be deduced by taking the limit fc —)• 0 
and s —>■ 0 which implies: 

d>(s) « 1 - STo-b ..., p(/c) « 1 - D|A:|“. (7) 


After proper rescaling of the waiting times and jumps 
the diffusion limit of the integral Eq.(|^ is obtained 
in the form of a space fractional Fokker-Planck equation 
(FFPE) 

^ Wo(a::, t\xo) = D-^—Woix, t\xo), (8) 

at o\x\°‘ 

with denoting the symmetric Riesz space fractional 
derivative which represents an integro-differential opera¬ 
tor defined as (33 : 


2 cos(7rQ!/2)r(2 — a) 

r -^^^dx' 

dx'^ J_^\x-x'\‘^-^ ’ 


(9) 


which has a particularly simple form in the Fourier space 


The first term on the RHS of the above renewal 
equation represents the survival of the probability mass 


F 


- gn 

_d\x\°^ 



-\k\<^F[f[x)]. 


( 10 ) 













3 


The total propagator of the process W{x,t\xQ) can then 
be obtained from Eq.([^. In the Laplace domain this 
equation has the form 


The Linnik PDF can be expressed in terms of elemen¬ 
tary functions only for a = 2, in which case it becomes a 
well-known Laplace distribution: 


W{x, s|xo) = Wo{x, s + r\xo) + -Wo{x, s + r|0).(ll) 

s 


Ps{x-,X,2) 



(17) 


In the case of Levy flights Wo{k,s\xo) 
Hence W{k,s\xo) is given by 


_D|fe|“-|-s ■ 


W{k,s\xo) 


gifcxo _j_ r 

-I- s -I- r 


( 12 ) 


and obeys the differential equation 

slT(/c, s|xo) - = 

= —D\k\°'W(k, slxo) — rW(k, slcco) -f (13) 

s 

The inverse transformationgives the FFPE describing the 
evolution of the total probability distribution: 

^^W{x,t\xo) = DJI^W{x,t\xo) - 

—rW{x,t\x{)) + r5{x), (14) 


with a zero mean and a variance Var[a;^] = 2A^. For 
a = 1 the closed-form expression for the corresponding 
Linnik PDF can be obtained (cf. Appendix]^ in terms of 

X oo 

special functions Si(a:) = / and Ci(a;) = - / 

0 X 

and in a scaled form reads: 

Aps(Ax; 1,1) = f ^ sin |a::|-Ci(|a;|) cos a;. 

\2 tt / 

(18) 

When passing to the analysis of the first arrival times in 
a subsequent Section, we note here that the result © 
has been obtained earlier in Ref. [23] for a discrete time 
counterpart of the resetting model. 


III. THE PROBLEM OF THE FIRST ARRIVAL 
TIME 


with initial condition lF(a;,0|a;o) = 5{x — Xq). Equation 
(141 is analogous to the Fokker-Planck equation defining 
a model of diffusion with stochastic resetting |5]. The 
difference lies in the fact that instead of a second or¬ 
der spatial derivative, characteristic of normal (Gaussian) 
diffusion, we are dealing now with a non-local fractional 
derivative, which describes Levy flights. Note that the 
model analyzed in this paper includes the other one as a 
special case, for a = 2. 

Having calculated the propagator in the Fourier- 
Laplace space, it is straightforward to obtain a character¬ 
istic function of the stationary distribution. For the sake 
of simplicity, we also introduce a length scale A“ = ^. 
By definition, the stationary PDF can be then derived 
from the relation 


Ps{k] A, a) = lim sW (fc, s|a;o) = 
s—>-0 

= rW'o(fc,s = r|0) = (15) 


The resulting function, Eq.(15l, is known as the Linnik 
distribution iiniiii], which is a special case of the family 
of geometric stable PDFs, approximating a distribution 
of normalized sums of i.i.d random variables 


N 

Sn — ^ ) Xi , 

i 


(16) 


where the number of terms N is sampled from a geomet¬ 
ric distribution, i.e. P{N = k) = (1 — Sum¬ 

mation of that type has been used, among others, in 
modeling energy release of earthquakes, water discharge 
over a dam during a flood, or avalanche dynamics |42j . 


For the stochastic process defined by Eqs.(Ip), a ques¬ 
tion of interest is in the estimation of the waiting time 
before the first event of a magnitude greater than a 
given threshold is observed. However, as has been dis¬ 
cussed elsewhere [44ti46] . the superdiffusive nature of 
Levy flights strongly influences the statistics of first pas¬ 
sage times over the threshold. In particular, due to long- 
range Levy jumps occurring with an appreciable proba¬ 
bility, the trajectory of the process may cross the thresh¬ 
old numerous times without actually hitting it. In conse¬ 
quence, the statistics of first arrival times at a predefined 
barrier is different from the statistics of first passages 
over it. 

Following Refs.fU [44], we introduce the first arrival 
time PDF pfa{t,x), which describes the distribution of 
times Tfa in terms of the integral equation for the prop¬ 
agator W{x, t|0, 0): 


t 

iy(x,^|0,0) = j dr pfa{T,x)W{x,t\x,T). (19) 
0 

The above formula can be easily interpreted: it simply 
states that the process which at time t finishes up at a:, 
has had to get to that point for the first time at some 
time T G (0,1). After that it could move freely until at 
time t, it came back to the very same point. The assump¬ 
tion of time-homogeneity {W (x, t\x, t) = W (x, t—T\x, 0)) 
explains a convolution operator on the RHS of Eq. ©. 
The function pfa is a probability density function of its 
first argument. The second argument denotes that the 
first arrival to a position x is evaluated. For readability, 
we skip D, r and a in the parameter list. From now 












4 


on, we also assume that the initial and reset positions 
coincide. 


By transforming Eq.(19l into the Laplace space a sim¬ 


ple algebraic relation is obtained: 

W{x,s\0) = Pfa{s,x)W{x,s\x). 


( 20 ) 


It is important to notice that W{x,s\x) ^ VL(0, s|0), as 
the resetting mechanism introduces space inhomogeneity. 
Our aim is to derive a formula for the mean first arrival 
time (MEAT) which can be obtained from pfa{s,x) as 
follows: 

(Tfaix)) = (21) 

We proceed by inserting the propagator, Eq. 0 , and 
the algebraic relation between the propagator and p/a(s)j 
Eq.(20l, into the formula for MEAT, Eq.([2l|). After 


straightforward algebraic manipulations we arrive at: 


/rr I fWo{x,s = r\x) 

= ? (ir„(»., = .| 0 ) - V = 

1 ^ Ps(0;A,a) _ ^ 


r \psix]X,a) 


( 22 ) 


Note that for simplicity we use a shortened notation 
Wo{x,t) = Wo(a;,t|0,0). Eq.(22) shows that the MEAT 
can be expressed either in terms of the Laplace transform 
of the propagator of the standard Levy a-stable process 
without resetting, or in terms of the stationary PDF of 
the process with the resetting mechanism switched on. 
This result is very general, since in the derivation no 
particular form of Wo{x,t\xo) has been assumed. 

We further focus on the special case of Levy flights. 
In general, the propagator of a Levy stable process can¬ 
not be expressed in terms of an elementary function of 
X. Representations in terms of the Fox functions m and 
in terms of the generalized hypergeometric functions |48j 
are known, but they are not useful in our case. We can, 
however, calculate Wo(a:, s = r\x) and deduce from its 
form the range of the stability parameter a that guaran¬ 
tees finiteness of the evaluated MEAT: 


A. Asymptotic behavior 

The average {Tfa{x)) cannot be expressed in terms of 
elementary functions for arbitrary a. Nevertheless, we 
can learn something about its behavior for large and 
small distances a: to a target. By taking a well-known 
expression for the asymptotic expansion of a-stable dis¬ 
tributions m and transforming it to the Laplace space, 
or otherwise, directly expanding: 


1 


D\k\°‘ + s 


g. {-D\k\°^y 

n—0 


on+l 


(24) 


and transforming this back from the Fourier space term 
by term, we obtain an asymptotic expansion of the prop¬ 
agator Wo in the Laplace space (see also [IHl EHl for more 
formal derivations): 


Wo(x,s|0,0) = 


1 ,D"r(na + l) 

-> (-1)"+^ sm(-na) ' . 

TT ''2 ^ gn+l^na+l 


(25) 

This expression is correct for a G (1,2). For a = 2 we 
don’t need the asymptotic expansion since in this case 
we have a closed-form expression: 

Wo(a;, s|0,0) = —(for a = 2). (26) 

2v Ds 


One can easily verify that Eqs.( 22|^ ) together with 
Eq.(|^ give the same result as the one derived in [H]. 
We truncate the series at the first term and so obtain the 
large x behavior of the MEAT: 


{Tfa{x)) CX 


rOl+l 


;1 < a <2 
; a = 2 


(27) 


We may also expand the MEAT around x = 0 using 
the known expansion of the Linnik distribution |49[ I50j . 
This leads to: 


(Tfaix)) 


a sin — 


1 




2sin^^^r(a)rsDi 


0{x 


2a-2 


)• 

(28) 


Wo{x,s\x) = ^ 


dk / = 


R 0 

r(i)r(i-i) 

iraDi s^~i 


1 


a sin -Do, 


(23) 


For any a; ^ 0, the propagator Wo(a:,r|0) is finite, since 
it is an integral of an oscillating function with an am¬ 
plitude decreasing to zero, and can be rewritten as an 
alternating series. We therefore conclude from Eq.(23) 


that the MEAT diverges for a < 1 and remains hnite for 
1 < a < 2. That apparent finiteness of the MEAT in case 
of Levy flights is rather surprising, taking into account 
the discontinuous character of superdiffusive trajectories 
and thus the possibility of overshooting (i.e. jumping 
over the target). 


IV. OPTIMIZATION OF THE MFAT 

Given a distance to a target x, one could be tempted 
to determine the optimal search kinetics of this location. 
We choose MFAT as an objective function, and mini¬ 
mize it in the space of parameters (r, a). We will de¬ 
note derived parameters of the efficient strategy as r*{x), 
a*{x), respectively and the corresponding optimal MFAT 
as T* (x). 


A. Fair comparison 

Since we want to compare Levy flights with different 
stability indices a, it is important to carefully choose the 
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parametrization of the family of jump distributions. One 
commonly used is ipik) = which in our case means 

fixing D = 1 for every a. Alas, this choice is very arbi¬ 
trary and based on simplicity of a characteristic function 
for symmetric stable distributions. As an alternative op¬ 
tion, we propose here a straightforward and consistent 
approach based on fractional moments. Let us define a 
random variable to be a position of the process with¬ 
out resetting at time t = 1 (this fixes the time unit). The 
p-th fractional moment may be expressed as 

Ag = (|ean = ^“/(a,p), (29) 

where the condition p < a has to be satisfied in order 
for the fractional moment to be hnite. Function f{a,p) 
is known and reads [51] 


/(a,p) 


2P+ir(£±i)r(-f) 

av^r(-i) 


(30) 


We want to keep Aq constant (e.g. A = 1) so our D will 
depend on a and p. The most natural choice of p in our 
case is p = 1 since it does not exclude any solution (in 
line with findings of Section III, we refer to cases with 
a > 1 assuring finiteness of MFAT) and it induces an Li 
norm that is commonly used in many applications. This 
choice leads to the expression: 


D{a) 



( 31 ) 


In the following we will refer to this method of compari¬ 
son, based on the choice p = 1, as the “fair comparison”. 
This is in contrast to the “naive comparison” based on 
the simplicity of characteristic function {D = 1). 


B. Asymptotic analysis 


From the asymptotic behavior of the MFAT several 
conclusions may be drawn: The prefactor in Eq.(28) is 


bounded for a € (1,2]. Consequently, for given non-zero 
r and D it is always possible to hnd x small enough, 
so that a = 2 minimizes the MFAT. In other words, 
Brownian motion is expected to be the optimal strategy 
at small distances to the target. In contrast, as it can be 
inferred from the asymptotic behavior, Eq.(27), for large 


enough distances the MFAT increases with x much faster 
for a = 2 than for a < 2. In this case, the Levy motion 
with a < 2 minimizes the MFAT, thus indicating a more 
efficient kinetics of space exploration to detect a target. 


C. Random distribution of target sites 

In many natural scenarios, living organisms navigate 
to unpredictable or randomly distributed resources. In 
other words, positions of the ’’target” is not precisely 


known. How is the kinetics of random search with re¬ 
setting affected by the location of targets in an unknown 
environment? In order to address this point, we further 
explore the MFAT under the constraint that the searcher 
knows only the mean (expected) distance to the target. 
Accordingly, instead of a fixed x in the evaluation of the 
MFAT, we use the PDF that satisfies the maximum en¬ 
tropy principle, i.e. a Laplace distribution p(x) of target 
positions is assumed. The MFAT in this more general 
setting can be calculated by averaging over possible dis¬ 
tances: 

OO 

(T/a(At)) = J dx{Tfa{x))p{x) = 

— OO 

OO 

= ^/ dx{Tfa{x))e~'^. (32) 

— OO 

Even though {Tfa{Xt)) is a different function from 
{Tfa{x)), for readability we keep the same symbol for 
the MFAT averaged over the distribution of targets and 
denote that by use of a different argument, only. 

As explained in the following example, such averaging 
over random distances to a target leads to modihcation of 
the MFAT and becomes crucial for the optimal strategy 
planning. Let us assume Brownian diffusion a = 2 with 
the Laplace PDF of target positions characterized by the 
mean distance to the target (lx]) = Aj. In that case the 
MFAT is given by the formula: 

iTfaiXt)) = l^, ( 33 ) 

' \ -L 


where A = • Clearly, the MFAT is finite for A > A* 


and optimization of Eq.(33l yields the value of the re¬ 
setting frequency r 2 (At) = If a searcher does not 

know the distribution of target locations, but was able 
to estimate via several measurements the mean distance 
to the target, (jxj) « At, he might be prompted to use 
that fixed position for further optimization of the MFAT, 


{Tfaix = At)). The derived optimal r*, see Eq.(Bl|, 


when applied to the system with Laplace distributed 
distance-to-target, would then lead to an infinite MFAT. 
This apparent inconsistency demonstrates that for the 
proper minimization of arrival times, the actual form of 
distance-to-target distribution p{x) is indispensable. 

It can be easily shown that for heavy-tailed distance- 
to-target distributions, the Brownian strategy always 
gives an infinite MFAT. In contrast, strategies with Levy- 
distributed jumps (a < 2) may provide efficient algo¬ 
rithms for searching, for which the MFAT remains finite 
as long as the pix) distribution is characterized by a fi¬ 
nite variance. A simple example illustrating this case 
is optimization of the MFAT given by Eq. (32) with the 


Student’s t-distribution of distances to the target: 


p{x) = 






1 + 


i^A? 


(34) 
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In this case the integral in Eq.(32) is convergent iff con¬ 
dition a < ^ — 1 holds. Numerical integration of Eq.(32) 
for V = 2.7 and i/ = A leads to MEAT functions displayed 
in Figs[^and[^ 


D. Scaling 

Optimal parameters r*{x), a*{x) and optimal MEAT 
T*{x) depend on x and D. For the sake of simplicity, 
from now on we fix D. It will be useful to take advantage 
of dimensional analysis to calculate the scaling behavior 
of the optimal r* and MEAT for a given a. Let r* (x) and 
T*(x) be the optimal r and the optimal MEAT for fixed 
X, a and D. Up to an arbitrary multiplicative constant, 
the only combination of x and D that has the dimension 
of time is t This leads to the following scaling 

equations: 


same r, at least in the vicinity of the optimal r*. More¬ 
over, in this limit Levy flights become more resilient to 
changes in r, especially in the range r > r*. 

Results displayed in Figj^ have been further analyzed 
to derive minimal values of the MEAT with respect to 
a pair of parameters (a, r) for different values of a dis¬ 
tance to a target, x. Consecutive Figand Figj^show 
outcomes of the optimization procedure described in Ap¬ 
pendix]^ for the cases of the immobile target located at 
a distance x, and the target with position described by 
Laplace distribution with an average distance to a target 
At, respectively. 

No qualitative difference in the derived optimal MEAT 
values has been found between the naive and the fair 
comparison. We therefore present results of the numer¬ 
ical optimization of {Tfa{x)) and {Tfa{Xt)) for the fair 
comparison only. 


r^(x) = r*(i)x° 

r*a{x) = 


(35) 


One easily verifies that these equations hold, by calcu¬ 
lating the derivative of the MEAT (Eq 22) with respect 
to r, comparing it to 0, and rewriting the corresponding 
equation such that it contains only a function of rx“. 
Scaling equations ( |3^ also imply similar relations to be 
fulfilled by 


As expected, for small x (A() Gaussian diffusive motion 
{a* = 2) is the optimal searching strategy. With growing 
distance to a target x (or At) the minimum of (Tfa) be¬ 
comes shallower, up to some point x* « 10.8 (Aj « 3.25), 
beyond which Gaussian diffusion is not efficient anymore 
and the optimal stability index switches to values a* < 2. 
Corresponding values of bifurcation points x* and Aj 
have been obtained by means of a numerical optimiza¬ 
tion procedure and are marked in Figs|3pl with a cross 
sign. 


T:iXt) = T:{l)Xf 


(36) 


The above relations are used in a numerical algorithm for 
optimization, as explained in details in the Appendix |B| 


E. Results 


A comparison between analytical prediction, Eq. (221, 


and numerical stochastic simulations has been performed 
and the results are displayed in Figj^ demonstrating a 
perfect agreement between both approaches. Addition¬ 
ally, Figj^ presents the analytically derived MEAT func¬ 
tions in 2-dim (a, r) parameter space. 

The MEAT diverges as r —0 and r —>■ oo (cf. Figs. 
HI- Accordingly, a minimum of the MEAT with respect 
to r can be found in the interval [0, oo) and its position 
depends on the stability index a characterizing underly¬ 
ing diffusive process. 

For small x MEAT values are systematically higher for 
non-Gaussian diffusion (a < 2) than for the Gaussian 
case and the same resetting rates. Also, as displayed in 
Fig|2 the MEAT has a more pronounced, deeper mini¬ 
mum in function of r for Levy diffusion with heavier tails 
(i.e. lower a’s), which suggests that the Gaussian strat¬ 
egy is more robust to variations of r. This is, however, 
no longer true for large x (cf. Figj^. In that case MEAT 
values for a = 2 are higher than that for a < 2 and the 


The described scenario of the continuous transition be¬ 
tween the Gaussian and non-Gaussian optimal strategies 
is qualitatively similar to the one investigated in Ref. [52j . 
In that article yet another variant of a one-dimensional 
Levy flight search strategy has been analyzed: The opti¬ 
mization of the random search for targets has been per¬ 
formed with respect to the average over inverse search 
times. The model has been enriched with a nonzero drift 
term (representative of an external bias or former expe¬ 
rience of the searcher) and no resetting mechanism has 
been included. Despite these differences, their plot of the 
optimal a* as a function of the initial position x (see Fig. 
3, Ref. [52]) at vanishing drift strength looks very similar 
to ours findings in Figs. 3 and 4: There exists a finite 
region (of relatively small x’s) in which the Brownian dif¬ 
fusion is the most efficient strategy and the optimal a* 
is always larger than 1. It seems that these observations 
are generic features of the analyzed optimal first arrival 
times. 

We have also investigated the impact of the heavy¬ 
tailed distribution of distances to the target on the effi¬ 
ciency of the searching. The analysis of the the optimal 
MEAT performed in this case is illustrated in Figsj^and 
1^ Presented plots indicate that the heavy-tailed distri¬ 
bution of distance-to-target excludes Gaussian diffusion, 
a = 2, from the set of possible optimal search strategies. 
Moreover, in line with the analysis of Section |IV C[ for 
Levy flights a condition a < v —1 has to be met in order 
to perform a successful search with a finite MEAT. 
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FIG. 1. Comparison between MFATs obtained by numeri¬ 
cal integration of the analytical formula (221 (lines) and by 
averaging over N = 10® realizations of a simulated process. 
Different lines (from the top to the bottom) correspond to 
a = (1.4,1.6,1.8, 2). For the sake of simulation not only time 
has to be discretized (St), but also a finite target size is needed. 
For each a the target size is chosen separately to match the 
analytical result at a; = 1, r = 1. The same target size is fur¬ 
ther used across different values of x and r. Estimated error 
bars are smaller than the markers used in the plots and hence 
have not been displayed. 


V. CONCLUSIONS 


Not only animal foraging patterns, but also memory 
retrievals of humans [53] and fluctuations of their sponta¬ 
neous activity [54] exhibit scaling statistics. The problem 
devised in this paper models mechanism of stochastic re¬ 
setting, or relaxation of a diffusive searching process to a 
predefined threshold, and as such can be well adapted to 
many natural scenarios of exploration processes such as, 
e.g. quests for food in a given territory |^, translocation 
and recruitment of repair proteins seeking for a disrupted 
DNA strand to be repaired [56] . optimal computer-aided 
web search [T2| or statistics of recall periods in retrospec¬ 
tive memory [ 55 ] . 

The efficiency of a search may be defined and analyzed 
by use of different measures, like e.g. the number of 
encounters of searchers and targets per unit of time, the 
mean inverse search time [52] . or the exploration range 
of space per unit of time. Here, we have focused on the 
efficiency measure expressed by the mean time to reach 
an immobile target, the MFAT. 

The first arrival time statistics has been analyzed for 
the one-dimensional problem with a constant resetting 
rate r. The acts of trajectory relocation have been as¬ 
sumed independent from the free (super)diffusive motion 
described by Levy jumps with the exponent 0 < a < 2. 
Despite the discontinuity of trajectories, typical for Levy 
flights, the MFAT remains finite iff a > 1 with a rich 


FIG. 2. The MFAT as a function of the parameters {a, r) 
for different values of the distance to an immobile target x = 
(1,10,20,100). Gontour plots beneath the surfaces help to 
guide an eye towards the minimum. 



X log 10 M 



FIG. 3. Optimal parameters {a* ,r*) and the MFAT as func¬ 
tions of the distance to a target x. 


characteristics of optimal (minimal) times T*{x). By use 
of the designed optimization method (Section IV), we 
have been able to derive the optimal parameters r*(x) 
and a*{x) for the range of target positions x. We have 
shown that the randomized distribution of targets with 
some average distance to a target results in a severe re¬ 
duction of distances for which Gaussian search remains 
the optimal strategy. Moreover, our analysis of optimal 
searching times for exponential distribution of distances 
to a target (Section [IV Cj ) clearly indicates that not only 
first moment of that distribution but rather its actual 
form is needed for a proper optimization planning: an 
optimization procedure based solely on the information 
about the average distance to a target would result in the 
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FIG. 4. Optimal parameters (a*,r*) and the MFAT as func¬ 
tions of At. The target position is a random variable with a 
Laplace PDF of distances and an average distance-to-target 
At. 


FIG. 6. Optimal parameters {a*,r*) and the MFAT as 
functions of At. Target position is given by a Student’s t- 
distribution, Eq.(34|, with v = 2.7 and an average distance- 
to-target At. 


VI. ACKNOWLEDGMENTS 




log 10 (A,) 



FIG. 5. Optimal parameters {a*,r*) and the MFAT as func¬ 
tions of At. The target position is given by a Student’s t- 
distribution, Eq.(34|, with = 4 and an average distance-to- 
target At. 


optimal r* leading to an infinite MFAT. 

Altogether, the proposed optimization scheme and 
scaling analysis can be further exploited, e.g. for two- and 
three-dimensional searching scenarios. Another plausi¬ 
ble modification of the proposed procedure could be an 
implementation of Levy walks, with coupled space-time 
distributions, or truncated Levy flights, penalizing very 
long jumps. 
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Appendix A: Linnik distribution 


The derivation of the Linnik PDF, expressed in terms 
of ps(a:. A, 1) in Eq.(18) proceeds as follows. Since 
Ps(a;, A, 1) as a function of x is even, without loss of gen¬ 
erality, we can assume that a: > 0. 


f{x) =TTps{x,l,l) = \J 

K 


dk 


^—ikx 


l + \k\ 


oo oo 


= J dk J ds cos {kx)e = 


= ds 


se 


= / dt 


te 


1+t^ 


(j r p ^21 J 


One can verify that g{x) is a solution of the equation 


g''{x) ^ g{x) = -, (A2) 

X 

which is a second order inhomogeneous linear differential 
equation with constant coefficients. We can easily solve 
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it by using the method of variation of parameters. Two 
constants in the general solution are calculated from the 
boundary conditions g{0) = f and lim g{x) = 0. The 

^ x—^oo 

solution reads: 




cos X + Ci(a;) sinx, 


(A3) 


which, after differentiation, leads to formula (181. 


Appendix B: Numerical scheme 


The optimization problem at hand could not be solved 
analytically. We have thus solved it numerically. Scaling 
formulas Eq. (35) allow for very fast numerical optimiza¬ 


tion, by reducing numerical calculation of the MEAT to 
one value of x for each a and r. The algorithm then 
proceeds as follows: For each a we perform numerical 
integration by use of the reverse Fourier transform of 
the Linnik distribution, Eq.(15), for a given value x, e.g. 


X = 1, and a few values of r. We fit a quadratic function 
to the calculated points and find the minimum of that 
function. Next we refine the interval of r values, center¬ 
ing it at the estimated minimum and, consecutively, we 


reduce its length. This procedure is repeated until the 
desired accuracy is achieved. We end up with a quadratic 
function which, by means of its vertex coordinates, de¬ 
fines our T*(l) and x*(l). Scaling equations, Eq.(35l, 
allow us to extend these results to arbitrary x. 


When we start the calculation for a new value of a, we 
face the problem of choosing a proper interval of values 
of r. Since we fit a quadratic function, it is important 
that the interval contains the optimal r. For this reason, 
we can make use of the optimal that was calculated 
in a previous step, for a value of a close to the new one. 
Accordingly, we choose an interval of r which contains 
’’’prev The formula for the optimal resetting frequency 
for the Brownian motion case is known [8] and reads: 





(Bl) 


with z « 1.5936. Therefore, when performing numerical 
analysis, we have started our calculations from a = 2. In 
the very last step we find, for each x, an a-parameter for 
which the smallest optimal MFAT is obtained. This is 
our global minimum. The numerical scheme used for the 
optimization of {Tfa{Xt)) is analogous. 
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